Zener tunneling isospin Hall effect in HgTe quantum wells and graphene multilayers 
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A Zener diode is a paradigmatic device in semiconductor-based electronics that consists of a 
p-n junction where an external electric field induces a switching behavior in the current-voltage 
characteristics. We study Zener tunneling in HgTe quantum wells and graphene multilayers. We 
find that the tunneling transition probability depends asymmetrically on the parallel momentum 
of the carriers to the barrier. In HgTe quantum wells the asymmetry is the opposite for each spin, 
whereas for graphene multilayers it is the opposite for each valley degree of freedom. In both cases, 
a spin/valley current flowing in the perpendicular direction to the applied held is produced. We 
relate the origin of this Zener tunneling spin/valley Hall effect to the Berry phase acquired by the 
carriers when they are adiabatically reflected from the gapped region. 



I. 



INTRODUCTION 



A large class of semiconductor devices is based on 
quantum mechanical tunneling of carriers through poten- 
tial barriers. This is the case of the Zener diode, which 
consists of a p-n junction where a strong enough elec- 
tric field induces interband transitions from the valence 
band of the p-type material to the conduction band of 
the n-doped material, see Fig. Ufa). The tunneling am- 
plitude is highly non linear in the applied field, and the 
tunneling current shows a breakdown-type behavior in 
the current-voltage {I-V) characteristics. The nonlinear- 
ity of the Zener tunneling makes this device very useful 
for semiconductor-based electronical. Interband tunnel- 
ing has been studied extensively in parabolic band-gap 
semiconductors, it is a paradigmatic example of non- 
adiabatic transitions, and it is known as the Landau- 
Zener tunneling^-. The most used model for studying 
the interband tunneling in parabolic semiconductors is a 
two-level system described by a Dirac-like Hamiltonian 
with a mass terrrPE!, see Fig. fl[b) . In this kind of ma- 
terials, the spin of the carriers typically plays no role. 

In this paper we are interested in analyzing Zener tun- 
neling physics in systems in which there is a correlation 
between the carrier's spin (or an equivalent degree of free- 
dom) and its direction of motion, i.e., systems is which 
chirality plays a role. In particular, we analyze two types 
of materials, HgTe quantum wells and carbon-based pla- 
nar heterostructures. These materials have in common 
that they can be described by 2 x 2 Hamiltonians and, 
therefore, it is possible to map the tunneling problem to 
the evolution of a two-level systenP. 

In the first case, HgTe quantum wells, we find that 
the Zener tunneling depends asymmetrically on the par- 
allel momentum of the carriers to the barrier, and this 
asymmetry is the opposite for each spin. We call this 
phenomenon Zener tunneling spin Hall effect. In these 
quantum wells the central region is an inverted band- 
gap semiconductor, such as HgTe, whose intrinsic strong 
spin-orbit coupling induces an inversion of the normal 
band progression of typical semiconductors, like the one 
used for the barrier material (e.g. CdTe). This kind of 
materials have come to the spot-light recently because. 
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FIG. 1: (Color online) (a) Schematic representation of a tun- 
neling process in a Zener diode in the uniform electric fleld 
approximation. T(ky) represents the transition probability of 
a quasiparticle with momentum ky from the p-doped valence 
band to the n-doped conduction band, (b) Wavevector-energy 
scheme for the two-band model. The diabatic (adiabatic) en- 
ergy levels are plotted in dashed red (solid blue) lines. Near 
the anticrossing region a diabatic e+ state can tunnel to the 
diabatic state e_. 



depending on the width of the central region, the sys- 
tem can undergo a quantum phase transition and become 
a topological insulatoiI^Hll. A topological insulator is a 
novel quantum state of matter that has metallic surface 
states inside the bulk energy gapPStlS. 

In the second case, graphene multilayers, we find that 
the Zener tunneling is also asymmetric with respect to 
the parallel momentum (except for monolayer), but the 
asymmetry changes for each carrier's valley index (in- 
stead of real spin). We call this phenomenon Zener 
tunneling valley Hall effect. Graphitic systems are also 
of great interest in condensed matter physics since it 
became possible to isolate monolayers, bilayers and in 
general multilayers of graphen^l-^*^^. p-n junctions of 
graphene have been created by gating locally these lay- 
ers and the transport properties of these heterostructures 
have been studied theoretically and experimentallj^^ ^■^. 
In particular, in bilayer graphene it is possible to open 
a gap in the spectrum by applying a voltage difference 
between the layers, and Zener tunneling is expected to oc- 
cur. Actually, it has recently been predicted that the I-V 
characteristics in bilayer graphene p-n junctions present. 



on top of the nonlinear Zener signal, some N-shaped 
branches with negative differential conductivity's. 

In both type of materials, the low-energy Hamiltonian 
can be expressed in terms of a pseudospin vector that 
multiplies the vector of Pauli matrices. In a tunneling 
process, the pseudospin vector undergoes a certain trajec- 
tory in the Bloch sphere and the carrier's wave function 
may acquire a Berry phase. We relate the Zener tran- 
sition asymmetry with the spin/valley-dependent Berry 
phase that the carriers acquire when they are adiabati- 
cally reflected from the gapped region. 

The paper is organized in the following way: In Section 
II we define the Hamiltonians that govern the properties 
of HgTe quantum wells and graphene multilayers. In Sec- 
tion III we map the tunneling problem to the time evo- 
lution of a two-level system and show numerical results 
for the different Hamiltonians. In Section IV we derive 
analytical expressions for the tunneling transition in the 
sudden and adiabatic approximations. In Section V the 
asymmetry of the tunneling amplitude as function of the 
momentum parallel to the barrier is explained in terms of 
the Berry phases that the carriers acquire upon reflection 
from the barrier. In Section VI we show the I-V char- 
acteristic curves for HgTe quantum well and multilayer 
graphene Zener diodes. We finish the paper in Section 
VII with a summary of our results. 



II. HAMILTONIANS 

A. HgTe quantum wells 

We study a HgTe quantum well confined by CdTe bar- 
riers. In bulk, and due to the strong spin-orbit coupling, 
HgTe is a zero gap semiconductor. When confined, HgTe 
is a normal band insulator for well thickness narrower 
than 63A and becomes a topological insulator for larger 
width^. For HgTe quantum wells grown in the (100) di- 
rection, the z-component of the spin is conserved. Near 
the band gap there are four relevant bands: the El bands 
that consist of the two spin states of the s-orbital, and 
the two spin states of the HHl bands which are a linear 
combination of p^ and pj^ orbitals. The low-energy effec- 
tive Hamiltonian for the two spin orientations, s^ = ±1, 
reads 



Hs, (k) = e{k)l + Mik)a, + A{kycry 



Szf^X^X } 7 



(1) 



where \i.= {kx, ky) is the in-plane wavevector of the carri- 
ers, k = |k|, e{k)=C-D{kl + fc2), M{k) = Mo + B{kl + 
ky), ax, (Jy, and Cz are the Pauli matrices and / is the 
identity. Eg = 2Mq is the band gap and A-D are parame- 
ters fitting the HgTe quantum wells^^. The product MqB 
determines the character of the insulator. For AIqB > 
the system is a normal insulator, whereas for MqB < 
a band inversion occurs and the system becomes a topo- 
logical insulator. The p-n structure of the Zener diode is 
described by adding to the Hamiltonian the appropriate 
scalar external potential of the form V{x)I. 



B. Multilayer Graphene 

In graphene and its multilayers, the low energy prop- 
erties occur near two non-equivalent valleys K and K' 
and the motion of the carriers depends on the valley 
where they reside. The role that the spin plays in 
HgTe quantum wells is played here by the valley index, 
r^ ± 1. Recently, it has been predicted that spin-orbit 
coupling in graphene opens a gap a nd the system could 
become a topological insulatop^SEni. However, this gap 
is very small and the occurrence of the quantized spin 
Hall effect would be observed at extremely low tempera- 
tures and in extremely clean samplea^D. Thus, neglecting 
spin-orbit coupling, the low energy properties of A^-layer 

ABC-stacke d mul tilayers are described, in general, by the 
Hamiltoniaii22l2ll 



H 



^ - / ,1. [cos{N(j)p)ax + sm{N(j)p)ay] + M^a, 
(-71) 



where p^ 



rirZ'r 



(2) 

= TzPy/p, 

In the previous expression the Pauli 



Here the notation is cos^p = Px/p and sint 



matrices act on the external layers for A^ > 1 and on 
atoms A and B of the unit cell in monolayer graphene. 
vf ^ Ix lO^ms"^ is the velocity of the carriers in mono- 
layer graphene22l and 71 '^ 0.3eV is the strongest direct 
interlayer hopping22. The last term in Eq. (l2|) opens a 
gap in the spectrum. In multilayer graphene this term 
represents an externally controlled potential shift in the 
chemical potential between the external layers. In mono- 
layer graphene, though, it is not possible to open a gap 
experimentally, but we are going to study this possibil- 
ity for the sake of completeness. Note, however, that at 
the surface of a 3D topological insulator there will exit 
a Dirac-like electron system'^'^ ^^ that, doped with mag- 
netic impurities, will develop a gap. The band structure 
of this surface state is governed by the same 2x2 Hamil- 
tonian than monolayer graphene, but with the cr-matrices 
referring to the real electron spin. 

As in the HgTe case, we describe the p-n structure 
of the Zener diode adding a scalar term, y(x)J, to the 
Hamiltonian of Eq. (pi) . 



HI. CONSTANT FIELD AND THE TWO-LEVEL 

SYSTEM 

In this work we describe Zener tunneling in the uni- 
form electric field model, V{x) — —Fx, see Fig. flja). In 
this approximation, it is possible to get the transmission 
across the p-n junction by mapping the problem into the 
evolution of a two-level system"* ^^. The key is that for an 
uniform electric field, F , applied in the a;-direction, the 
problem can be simplified if we use the momentum rep- 
resentation, X = idk^. With it, the Schrodinger equation 
corresponding to Eqs. (fTl) and ^ becomes 



zF^ = [Hs{k) ~ El]i;s, 



(3) 



where E is the energy and the index s stands for the spin 
or the valley index, depending on the system at hand. 
For each index, s — ±1, this equation is identical to 
the Bloch equation describing the dynamics of a spin- 
1/2 particle in the presence of a magnetic field, with the 
wavevector in the x-direction playing the role of time. In 
the uniform electric field model, the term [e(fc) — -E]/ for 
HgTe quantum wells in Eq. (l3| , or equivalently the term 
EI for multilayers, does not contribute to the interband 
transition and we drop it. 

In the limit k^ — )■ ±00 the eigenvalues of Hg are also 
eigenstates of a^, being v = x for multilayer graphene 
and V = z for HgTe quantum wells. Starting at k^ — > 
—00 from the low energy eigenvector (with eigenvalue 
a^, = —1) and tuning k^ from —00 to +00, the two level 
system traverses a level anticrossing^. Valence to conduc- 
tion interband transitions are described by the process 
in which a state that at k^ — —00 has negative energy 
evolves into a state that at k^ — + cx) has positive en- 
ergy. We have solved numerically Eq. ([3]) in an interval 

^x,min ^ ^x ^ '^x,max SUCU tnat, at f^x^max and tZx^mini 

the eigenvalues of (t^ are ±1. From the evolution of Eq. 
(|3|) we obtain the wavefunction at k^^max and from the 
square of its projection on the state with a^ = 1 we ob- 
tain the interband transition probability. 

In Fig. [2] we plot the Zener transition probability 
T(ky,s = 1,-F) as function of the wavevector ky of the 
incident particle for a HgTe quantum well and for mono- 
layer, bilayer and trilayer graphene. We plot the tran- 
sition probability for several values of the electric field, 
F=2, 5, 10 and 20i^o, with Fq = Mo/dM- Here Fq and 
dM are the electric field and length characteristic scales 
set by the gap of the insulator: d^.j = A/Mq for HgTe 
quantum wells and dM—fivpi^i — tt^t)^'^ for multilayer 
graphene. 

The symmetry of the Hamiltonian dictates that 
T(ky, s, F)=T(—ky, —s, F). Except for monolayer 

graphene, the transition probability has a maximum at 
a finite value of ky that depends on the sign of s. As a 
result, carriers with positive spin/isospin are mainly de- 
flected towards one y-direction when tunneling, whereas 
those with negative spin/isospin are deflected in the op- 
posite direction. The overall transition probability in- 
creases with the applied electric field, see Fig. [2] At 
small fields the spatial extension of the forbidden region 
becomes very large and the Zener tunneling amplitude 
goes to zero abruptly when F — ^ 0. This is the origin of 
the switching behavior of the Zener diodes. For moder- 
ate applied electric fields the asymmetry in the angle of 
incidence also increases with the field. 



IV. SUDDEN AND ADIABATIC 
APPROXIMATIONS 
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FIG. 2: (Color online) Zener tunneling in the constant field 
model as a function of the parallel momentum of the incident 
particle for HgTe quantum wells and multilayer graphene. 
The results correspond to spin/valley s—1 (for s = —1 equiva- 
lent results are obtained, but specularly reflected with respect 
to ky = 0) . The direction of the arrows indicate the evolution 
of the curves when increasing the electric field. 



calculations expand the solution of the Hamiltonian in a 
diabatic or in an adiabatic basis, see Fig. [lib). The first 
case is suitable for an unperturbed Hamiltonian that can 
be diagonalized in a diabatic basis, where the carriers 
evolve with probability one from the valence to the con- 
duction band. We then calculate the first correction to 
perfect transmission in the sudden approximation, treat- 
ing the rest of the Hamiltonian in first order perturbation 
theory. When the Hamiltonian is such that the tunnel- 
ing transmission form valence to conduction band is very 
small, it is more convenient to use the adiabatic basis as 
the unperturbed one. In the adiabatic basis the carriers 
are completely refiected at the barrier, and the tunneling 
probability can be obtained as first order perturbation to 
the adiabatic Hamiltonian. 

The Hamiltonians of the systems we are studying can 
be written in the form 



H = e(k)h(k)-a. 



(4) 



In order to shed some light on our nunrerical results, we 
have solved Eq. (Is]) analytically in the limit of small par- 
allel momentum ky and large electric field. The analytical 



This equation defines a wavevector-dependent unitary 
pseudospin vector /i(k). From the form of this Hamil- 
tonian, the expectation value of the vector of a matrices 



is either parallel, in the conduction band, or antiparallel, 
in the valence band, to the pseudospin. In the absence 
of gap, Mq = 0, carriers approaching perpendicularly to 
the barrier, ky — 0, should conserve the pseudospin. In 
agreement with the Klein paradox, that implies perfect 
transmission for gapless monolayer and trilayer graphene 
and perfect reflection for gapless bilayer graphene^^. This 
conservation of the pseudospin at fcj, = for massless 
Hamiltonians would help us to choose a diabatic or adi- 
abatic basis as the starting point in perturbation theory. 



A. HgTe quantum wells, diabatic basis and sudden 
approximation 

In HgTe quantum wells, the pseudospin has the form 

J^HgTe ^ (^^^^ ^^^^^^ Mo+ Bk^) / y^ {Mo+ Bk^)^ + A^k^ . 

(5) 
For Mq = and B = 0, the pseudospin takes the form: 
{kx,kySz,0)/\k\, and the eigenfunctions of the Hamilto- 
nian of Eq. ([l]) are chiral. In this limit, the Klein para- 
dox dictates ^^ that the tunneling amplitude at ky = is 
unity. For finite values of Mq and B, the Klein paradox 
does not apply exactly, but the transmission probability 
at large electric fields and small values of Mq and ky is 
close to unity. Therefore, it is convenient to work in the 
diabatic basis. 

In the natural units of the problem, x = k^dM, y = 
kydM, £ = F/Fq and B = B/{Mod\j), the Hamiltonian 
is written as: 



ox 



iHo + V)ip ,with 



Bx^ 



-Bx- 



V-- 



hBy^ 
iy 



-ly 
-l-By\ 



The eigenvectors of i/o define the diabatic basis, with 



eigenvalues e±{x) = ±^\/\ -\- B'^x'^, see Fig. [ub). We 
represent the corresponding wavefunction in the diabatic 
basis as 



V'(x) = C_(x)e-^"'(^)|- > +C+(a;)e+'"^(^)| + 



>, 



(7) 



where Hq\± >= £±|± > and uj{x) = / £j^{x')dx' . We 
define the transmission t and reflection amplitude r as 
follows: assuming that C^{~oo)— and C+(— oo)= 1, 
then r ~ C_(cx)) and t — C+(oo). Plugging the wave- 
function, Eq. (It]), into the Hamiltonian of Ec^. (l6|, we 
get 



dxC± — ±—ToC± 



T± Czr , where 



(8) 



ro=(i + Bj/2) 



Bx^ 



y x"^ + B'^a 



and 



— IIIp=F2i0(a;) 



r±=Te 



1 + By'^ 



V2 



B^3 



+iyH£s^ 



B/2 
l + Bx^ 



(9) 



The amplitudes r and t are obtained from the asymp- 
totic solution of Eq. ([8]) with the appropriate boundary 
conditions. For large values of £ it is possible to get an 
analytical expression for the transition. To lowest order 
in l/£, the asymptotic form of C_ is obtained by substi- 
tuting C+{x) = 1 in Eq. ^, 



C_(c5o) = - 



+CX3 



T-dx. 



(10) 



For small values of B^£ this integral can be evaluated 
using the steepest descent method. From it we obtain 
the interband transition probability 



T{y,s^ 



£ 



1 



i+Sy 



B3£2 



£B 



.s.(l- 



2' 



.(11) 



Recovering previous units, the maximum of the transition 
probability to lowest order in FB'^ occurs at a wavevector 



Ky - 



FB (^ 5 B , ^ 



(12) 



(13) 



being the maximum transition probability 

J^max - 1 PA' 

In Fig. [3] we compare the transmission at ky=0 obtained 
numerically from Eq. (|3| with the one obtained from 
Eq. ( 11 ). The quality of the approximation is good, spe- 
cially for strong electric fields. Equations ( 12 ) and ( 13 ) 



explain qualitatively much of the results presented in the 
first panel of Fig. [2] i) for each spin orientation the tran- 
■ (6) sition probability is asymmetric with respect to ky, ii) the 
asymmetry increases with the field, iii) the sign of ky^ de- 
pends on the product s^B (note that when B — there 
is no spin Hall effect), iv) the asymmetry is present either 
for MqB > or MqB < 0, i.e., irrespective of whether 
the quantum well is in the trivial or in the topological 
phase, and v) the overall transition increases with the 
electric field. Moreover, Eqs. (12|-(13) describe quanti- 
tatively the results in the case of large F. For example. 



for F/Fo^lO we get T„a^=0.69 and k^ = 0.45, results 
that are in rather good agreement with the numerical 
ones presented in Fig. [2j 

A spin dependent transmission has been also predicted 
to occur at the interface between a HgTe quantum well 
and a metaP^l. In this case the asymmetry is related to 
localized states at the interface. 



B. Monolayer graphene, diabatic basis and sudden 
approximation 

The pseudospin vector for graphene has the form 

/J™ = (vFPx. VFT,py,Mo)/jM§+vj,p-', (14) 
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FIG. 3: (Color online) Comparison of the transmission prob- 
ability at ky — obtained by solving numerically Eq. ([3| 
with the analytical result obtained in first-order perturba- 
tion theory: Eq. (Ill for HgTe quantum wells, Eq 
for graphene, Eq.(28l for bilayer graphene and Eq 
trilayer graphene. 



and for gapless monolayer graphene the Klein paradox 
applies exactly. Therefore, in order to study the tunnel- 
ing when Mq ^ 0, it is convenient to work in the dia- 
batic basis and, using natural units as before, we write 
the Bloch-likc equation as 



ox 



X 
X 



1 

iyTz 



-WTz 

-I 



^. (15) 



Note that in the reduced units, hvp plays the same role 
as A. The eigenvalues of the first term of Eq. ( 15 1 define 



the diabatic basis. Using this basis, the wavefunction 
takes the form 



^{x) = C^{x) 



1 



1 
-1 



e'^ +C+ix) 



1 



and the coefiicients C_ and C+ satisfy 



d..C^ = -^{l-iyT,)e'^C±. 



(16) 



(17) 



The reflection amplitude in the sudden approximation, 
valid in the f — >■ oo limit, is then 



r(l-iyT"z) 



'2edx, 



(18) 



being the transition probability, in the original units, 

T{ky) = 1 - ^(M2 + h^klvl) . (19) 

In agreement with the numerical results, we get that the 
transition probability is symmetric in ky and indepen- 
dent on the isospin r^ (see Fig. [2]). Note that, for mono- 
layer graghene, the transition probability can be obtained 
exact l-y 



TT91 



T{ky) 



(20) 



and Eq. (19) corresponds to the first term in the 1/F 
expansion. Fig. [3j second panel, illustrates the quality 
of the sudden approximation at large values of the electric 
field. 



C. Bilayer graphene and the adiabatic 
approximation 

For bilayer graphene, the pseudospin has the form 



"fo , T^.fA I . I i\,r2 _L !iF:^4 



-{pi - Py), —2p,PyT,,Mo /WM2 + -^p- 

.71 71 / V ^1 

(21) 
Gapless bilayer graphene, Mq = 0, is also chiral, and the 
pseudospin is (fc^ — ky,2kxkyTz,0)/k'^. Holes impinging 
perpendicularly to the barrier from the left have opposite 
pseudospin than electrons moving to the right, and the 
same pseudospin than holes refiecting from the barrier. 
Therefore, the tunneling probability for Mq = and ky — 
is null. When Mq ^ 0, the transition probability at 
fcj, = is still small, see Fig. [2] and it is thus more 
appropriate to work in the adiabatic basis. 

The Bloch equation for bilayer graphene has the form 



ox 



where now £ 



1 

{x + iyTzf 



[x ~ iyTz 
-1 



hvf 



(22) 



Mo VMrn' 
The Hamiltonian Hb in Eq. ([22|) defines the adiabatic 



basis with eigenvalues e{x,y) = ±-\/l + [x"^ +y'^Y and 
eigenfunctions 



7/i_ 



- sin (61/2) 
cos (61/2) e"^ 



and ip^ 



cos (6'/2) 
sin (6'/2)e^'^ 

■ (23) 
In order 






where cos [9) — I/|e| and (j) = tan ^ ( 

to solve Eq. ( |22| ) we consider the general solution 



Ip — aie 



"-0- +026* 



-^, 



+ j 



(24) 



with uj{x,y) = /„ e{x' ,y)dx' . The coefficients ai and a2 
satisfy 



Oi 



—a2T{x, y) — iai cos 



2jdx 



d_ 

dx 

-—a-z = -aiT*{x,y)-ia2sin'^ i o ; a ' 
ox \2 I ox 



(25) 



with 



(jpvl iyjix"^ - y^) + x{x'^ - ?,y'^)/e{x, y) 
{x"^ +y^)e{x,y) 



(26) 



Now we take the adiabatic hmit, i.e., we consider that 
the probabihty to undergo a transition from one adia- 
batic state to another is neghgible, 02 ^ 1. Then the 
transmission amphtude is 



tiy) 



T{x,y)dx. 



In the hmit y — > this integral can be evaluated follow- 
ing the methods presented in Refs. j6|38j and we obtain 
(back in physical units), see appendix A, 



4^^ 



A£o yj^oTl 



T{ky,T,) w ^e '"^"^" "^"'f' sin^ici 



Mo VMmA 



F hvp J 



X 1 - C2kyTz 



1/3 



(hvp) 



4/3 



Mo J (Mo7i)2/3 



(28) 



where ci ~ 1.23 and C2 ~ 3.19 are numerical factors. 
Some comments on Eq. ( 28 ) are in order: i) as shown 



probability at finite values of F, A/q and ky. These ze- 
ros appear because the bilayer graphene Hamiltonian is 
quadratic in the momentum and, for each energy in the 
gap region, there are two decaying states that interfere 
under the tunneling barriep22l. 



We write the Bloch equation as the sum of a diabatic 
term plus a perturbation. 



.d 1 

ox t 


'{0 x^\,{ 1 {x~iT^y)^-x^ 
\x^ ) ^ \{x+iT,y)^-x^ -1 



and in the case of the trilayer graphene we have 

Mo (7?Mo)i/3 



£ = 



F 



hvp 



(30) 



(31) 



(27) The eigenfunctions of the first term of Eq. (30 1 define 



the diabatic basis. In this basis the wavefunction takes 
the form 



V'(x) 



C-{x) 



1 



C+{x) 
V2 



(32) 



V2 V-1 

and the coefficients C_ and C+ satisfy 
d,C^^T^3xy^C±+^{l±iyT,{y^-3x^))e-'^C^ . (33) 
To lowest order in l/£ and y, the reflection amplitude is 



3/4 



in Fig. |3J there is a reasonable agreement between the 
numerical results and the one obtained in the adiabatic 
approximation, ii) the transition probability is not sym- 
metric with respect to ky, but it is so with respect to 
the product kyT^, iii) the maximum transition occurs at 
finite ky, iv) there is an oscillatory term in the trans- 
mission amplitude that produces zeros in the tunneling T{ky,Tz) — 1- 



r = C_(+oo) ~ -- / (1 -i3x^2;T2)e*"2 £ dx 



lV'Ml/4)„,- 



e 8 



1/4 

3r^2/ 



r(3/4)_,3. 



(34) 



^£ ) 23/4 " y^j " '-^ 21/4 

and the transmission probability, in the physical units, is 

M271 r2(i/4) Mofc^ r(i/4)r(3/4) 



{hvpFfl'^ 23/2 



%/2 



(35) 
This expression agrees remarkably well with the nu- 
merical results for large electric fields, see Fig. |3] Also, 



Eq. (35|) explains qualitatively the dependence of the 
tunnel: 



ing probability on the wavevector ky. 



D. Trilayer graphene, diabatic basis and sudden 
approximation 



BERRY PHASE AND LACK OF 
REFLECTION SYMMETRY 



In ABC-stacked trilayer graphene the pseudospin uni- 
tary vector is 



(f(4p? 



3p,p2)^^(3p2p^ 



4p3)r„Afo 



Ml -f ^p6 



(29) 



For massless trilayer graphene the pseudospin reduces to 
(4fc3 _ Skxk'^ , {3kyk^ - 4fc3)r^,0)/fc3 and the eigenvec- 
tors are again chiral. Because the pseudospin rotates 6tt 
when the wavevector rotates 2tt around k = 0, the tran- 
sition probability at ky—0 is unity for massless trilayer 
graphene. In Fig. |2]we see that, even for Mq y^ 0, in the 
limit of large electric field the transition probability at 
small ky is near one. Therefore, it is appropriate to work 
in the diabatic basis and use the sudden approximation. 



The question that remains is the physical origin of the 
asymmetry, for a fixed spin/valley, of the tunneling am- 
plitude as a function of ky. The asymmetry is not re- 
lated to Chern number associated with the chirality of 
the massless, Mq = 0, HamiltoniansSH. Although HgTe 
and monolayer graphene share the same Chern number, 
in monolayer graphene the transition amplitude is sym- 
metric with respect ky, whereas it is not so in HgTe quan- 
tum wells. 

We associate the asymmetry with the winding of the 
expectation value of the pseudospin /i(k) when a carrier 
is adiabatically reflected by the tunneling barrier. This 
is related to the sign of the Berry phase acquired by the 
carrier's pseudospin in this process. 

Consider a quasiparticle moving in the valence band 
in the presence of a constant electric field, V{x) = —Fx. 
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FIG. 5: (Color online) It-^.±-V characteristics of monolayer 
graphene for two values of the built-in potential. Units are 

h = ^VoW/dM X 10*^ and eVb = MoL/dM- The inset shows 
schematically the Zener diode. 



FIG. 4: (Color online) In-plane projection of the trajectory 
defined by the pseudospin h{kx, ky) in an adiabatic reflection 



process. Solid lines correspond to trajectories with sky 
and dotted lines to trajectories with sky < 0. 



> 



This quasihole coming from x = — oo and moving to- 
wards the right has a momentum fc^; < 0. Upon arriving 
into the gapped region, it is adiabatically reflected from 
it back to a; = — c» with momentum k^ > 0. In the pres- 
ence of the electric field, the momentum kx is not a good 
quantum number and it is not conserved. In a semiclas- 
sical/adiabatic approximation the momentum is defined 
by the relation 



e{ka;{x),ky) - Fx^ E 



(36) 



In the reflection process the pseudospin h{kx,ky) de- 
scribes a trajectory on the Bloch sphere of radius unity. 
When the trajectory closes a circuit F in the unit sphere, 
the wavefunction of the carrier acquires a Berry phase 
equal to half the solid angle defined by the surface en- 
closed by the circuit F. 

In Fig. |4]we plot, for the different Hamiltonians stud- 
ied in this paper, the in-plane projection of the trajec- 
tories defined by the pseudospin h when the carrier goes 
from X = — oo to the barrier and is reflected adiabatically 
back to a; = — oo. In the case of monolayer graphene such 
trajectory defines an open line, both for T^ky greater or 
smaller than zero. Thus, for monolayer graphene there 
are no closed paths in the adiabatic process and there is 
no Berry phase associated with the reflection. The situa- 
tion is different for HgTe quantum wells. In this case the 
trajectory defines a closed circuit and there is a Berry 
phase associated with the adiabatic refiection. The sign 
of the Berry phase depends on the direction in which 
the closed loop is traversed by the pseudospin. It turns 
out that it has opposite sign for opposite signs of ky or 



Sz- Therefore, the sign of the Berry phase depends on 
the sign of the product kyS^- For graphene multilayers, 
A^ > 1, the pseudospin trajectory in the adiabatic reflec- 
tion process always defines closed paths that have oppo- 
site orientation for opposite values of the product kyT^. 
The dependence of the Berry phase on the product fc^s, 
being s the spin or the valley index, breaks the reflection 
symmetry in each index s and explains the asymmetry of 
the transmission for a momentum ky at a fixed index s. 



VI. ZENER TUNNELING CURRENT 

Finally, we calculate the tunneling current flowing 
through the Zener diode. We consider a p-n junction as 
the one sketched in the inset of Fig. [5] U represents the 
built-in potential induced by doping or electrical gates, 
and Vsd is the source-drain potential difference. L is the 
junction length. Within the Landauer approximation, 
the tunneling current for index s moving in the positive 
y-direction has the form 



Is,+ = 



/oo 
-oo 



w 

2^ 



T[ky,s,-{U + Vsd)]dky, (37) 



where ue is the Fermi-Dirac distribution and W is the 
transverse length of the p-n interface. Although in the 
uniform field approximation the transition amplitude is 
energy independent, the limits of the integral in ky de- 
pend on energy through the relation E = s^k^ — 0,qy)- 
The current flowing in the negative y-direction, Is.~, is 
obtained performing the integral in ky from —qy to 0. As 
the transition is dominated by small values of ky, see Fig. 
[2] and Ref. 1271 we approximate Qy by oo in the calcula- 
tion of the currents. For zero temperature the current 
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FIG. 6: (Color online) I-t-,±-V characteristics of HgTe quan- 
tum wells, for two values of the built-in potential. Panels (a) 
and (b) correspond to current flowing in the positive and neg- 
ative y-direction, respectively. In panel (c) we plot the excess 
of current in the positive j/-direction for electrons with spin 
up. Units are 7o = ^VoW/dM x 10^ and eVb = MoL/cLm. 
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FIG. 7: (Color online) I^^=i ±-V characteristics of bilayer 
graphene, for two values of the built-in potential. Panels (a) 
and (b) correspond to current flowing in the positive and neg- 
ative y-direction, respectively. In panel (c) we plot the excess 
of current in the positive y-direction for electrons in the valley 
T^ = 1. Units are 7o = ^VoW/dM x 10^ and eVo = Mq. 
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FIG. 8: (Color online) It^=i,±-V characteristics of trilayer 
graphene, for two values of the built-in potential. Panels (a) 
and (b) correspond to current flowing in the positive and neg- 
ative y-direction, respectively. In panel (c) we plot the excess 
of current in the positive y-direction for electrons in the valley 
T^ = 1. Units are 7o = ^VoW/dM x lO"* and eVo = Mq. 



gets the form 



e2 W 

Is,± — '^-T^sdTT- 
h ZTT 



iboo 



T[k 



y, • 



L 



{U + Vsd)]dky. (38) 



This current verifies the symmetries Is^+=I-s.- and 

In Fig. [Slwe plot the Ir^.±i-V characteristics for mono- 
layer graphene and different values of U. The curves 
present the breakdown-type behavior characteristic of a 
Zener diode. In the case of monolayer graphene the trans- 
mission amplitude is symmetric with respect to the mo- 
mentum ky and the current is equal for positive and neg- 
ative y-direction. 

In Fig. [6]we plot the I^^±i-V curves for a HgTe quan- 
tum well and different values of U. For spin up there is an 
excess of current in the negative y-direction. This effect 
is the opposite for spin down electrons. These results in- 
dicate the existence of a spin current perpendicular to the 
Zener barrier. This Zener tunneling spin Hall effect is a 
consequence of the asymmetry in the transition curves of 
Fig. [2] From Fig. [6] we obtain that the Hall spin current 
can be as large as 30 per cent of the electrical current. 
Because a HgTe quantum well may be a two-dimensional 
topological insulator under certain conditions, there is 
an extra contribution to the current in this case com- 
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ing from the spin-polarized edge states developed in the 
barrier region. However, its magnitude is always of the 
order of one conductance quantum ^ e^ /h or less, since 
the electric field diminishes ilPS. On the other hand, the 
Zener tunneling spin Hall current is proportional to the 
transverse length W, see Eq. (38 1, and increases with 
the electric field. 

In Figs. [7] and |8] we plot the I-V characteristic curves 
for bilayer and trilayer based Zener diodes. In both 
cases there are some oscillations on top of the non-lineal, 
N-shaped I-V curves. These negative differential con- 
ductivities appear for positive and negative y-directions, 
and they have their origin in the interference between 
decaying states in the energy gap regiorp2. In trilayer 
graphene the negative differential conductivity is even 
stronger than in bilayer graphene. At large source-to- 
drain voltage, the asymmetry of the tunneling ampli- 
tude is reflected in an excess of current in the negative 
y-direction with respect to the positive y-direction. For 
multilayer graphene, this effect is the opposite depending 
on the valley r^ . These results indicate the existence of a 
valley current perpendicular to the Zener barrier that is 
a consequence of the asymmetry in the transition curves 
of Fig. [| 



VII. SUMMARY 

We have analyzed Zener diode physics in HgTe quan- 
tum wells and multilayer graphene. In the case of HgTe 
quantum wells we find that, after traversing the barrier, 
a Zener tunneling spin Hall current is developed to the 
right of the diode. In the case of bilayer and trilayer 
graphene the Zener diode generates a valley Hall cur- 
rent. This effect is absent for the monolayer graphene. 
The magnitude and polarization of the Hall currents in- 
crease with the applied electric field. 

The tunneling current is obtained from the transmis- 
sion probability that is computed numerically in the con- 
stant electric field approximation. The origin of the Hall 
currents is the asymmetry of the transmission probabil- 
ity in the momentum ky perpendicular to the tunneling 
barrier. We have developed an analytical approximation 
for the tunneling transmission at small ky that agrees 
rather well with the numerical results. The physical ori- 
gin of the Zener tunneling asymmetry on ky is related to 
the Berry phase that the carriers acquire when they are 
adiabatically reflected from the tunneling region. 

In the case of multilayer graphene the Zener tun- 
neling valley Hall effect could be used for valleytronic 
application ;] ' '^^ 1. In an appropriated geometry, the 
asymmetry in the Zener tunneling should enable to spa- 
tially separate the carriers of each valley^, which could 



be useful to manipulate the valley degree of freedom in 
bulk garphene. 

The Zener tunneling spin Hall effect we predict to oc- 
cur in HgTe quantum wells could be used for electrical 
manipulation of the spin currents. The spin currents in 
the Zener device should be stronger than those occurring 
in diffusive systems, and they could be detected in non 
local electrical measurements in H-shaped structures^. 
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Appendix A 



In this appendix we evaluate the integral of Eq. (27 1 



In the limit y — > 0, the transition probability takes the 
form 



dxe-^^^^f-^-fz^iL^). (Al) 



l + x^ 



Vl + x^ 



In terms of Lo{x^y = 0), the integral has simple poles in 
the complex plane at uji — Jq^ dxVT-j-a?, where Xi = 
_l_g±i7r/4^ Expanding the value of uj near Xi we find 



Ljix,y = 0)-u,^-x^^\x-x,f/^ (A2) 



and we rewrite 



x{uj) 



3iy 



We solve this integral by closing the path around the 
lower half of the complex plane. This path encloses the 
poles %/2/2(±l,— i) and their associated branches. The 
integral then yields 



t{v) 



27r 



e--iHin(ci^)(l-C2£'i/32/), 



(A4) 



with r. - 1 /TlXiM and m - ^^'^ 

witn ci - J ^ 2 r(774) ™^ "^2 - r(273) • 
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